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ABSTRACT 

We have created a general methodology for calculating the wavelength-dependent light curves of 
close-in extrasolar giant planets (EGPs) as they traverse their orbits. Focussing on the transiting EGPs 
HD189733b, TrES-1, and HD209458b, we calculate planet/star flux ratios during secondary eclipse 
and compare them with the Spitzer data points obtained so far in the mid-infrared. We introduce 
a simple parametrization for the redistribution of heat to the planet's nightside, derive constraints 
on this parameter (P ra ), and provide a general set of predictions for planet/star contrast ratios as 
a function of wavelength, model, and phase. Moreover, we calculate average dayside and nightside 
atmospheric temperature/pressure profiles for each transiting planet/P„ pair with which existing and 
anticipated Spitzer data can be used to probe the atmospheric thermal structure of severely irradiated 
EGPs. We find that the baseline models do a good job of fitting the current secondary eclipse dataset, 
but that the Spitzer error bars are not yet small enough to discriminate cleanly between all the various 
possibilities. 

Subject headings: stars: individual (TrES-1, HD209458, HD189733) — (stars:) planetary systems — 
planets and satellites: general 



1. INTRODUCTION 

Probing the atmospheres of extrasolar giant planets 
(EGPs) by measuring their spectra is the paramount 
means to determine their physical and chemical char- 
acter (Burrows 2005). Such direct measurements com- 
plement the kinematic and orbital information obtained 
through the radial-velocity (RV) technique by which the 
vast majority of the EGPs have to date been discov- 
ered and studied. However, an EGP's spectrum and 
phase-dependent light curve can in principle reveal or 
constrain the molecular and atomic compositions, atmo- 
spheric temperatures, cloud properties, albedos in the 
optical, and the degree to which the heat absorbed on 
the dayside is redistributed to the nightside before rera- 
diation. The advection of heat and material from the 
dayside by jet streams and zonal winds will alter the 
dayside atmospheric temperatures and non-equilibrium 
compositions (Menou et al. 2002; Cho et al. 2003; Burk- 
ert et al. 2005; Iro, Bezard, & Guillot 2005; Showman 
& Guillot 2002; Guillot & Showman 2002), and might 
measurably shift the light curve with respect to the or- 
bital ephemeris (Cooper & Showman 2005; Williams et 
al. 2006). Dynamic meteorology could also introduce 
zonal banding, as seen in Jupiter and Saturn, and tem- 
poral fluctuations, and does influence the rate with which 
heat is lost from the inner core (Burrows, Sudarsky, & 
Hubbard 2003; Burrows et al. 2004), and thereby the 
radius of the planet and its evolution. Moreover, such 
redistribution affects the near- and mid-infrared emis- 
sions from the night side, and as a result will affect the 
interpretation of nightside data when they become avail- 
able. 

The planet/star flux ratios of wide-separation EGPs 
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(> 0.2 AU) are quite low (10~ 4 to 10~ 14 ) and vary 
widely as a function of wavelength and orbital separation 
(with the concomitant non-monotonic variations in geo- 
metric and Bond albedos) (Sudarsky, Burrows, & Pinto 
2000; Burrows, Sudarsky, & Hubeny 2004; Sudarsky et 
al. 2005; Burrows 2005). Nevertheless, space-based coro- 
nagraphic techniques can be designed with inner working 
angles and contrast capabilities that will eventually im- 
age such planetary systems in the optical and mid-IR 
and distinguish planet from star (Trauger et al. 2000; 
Trauger, Hull, & Redding 2001). However, the close-in 
EGPs with orbital semi-major axes less than ^0.1 AU 
will not be imaged separately any time soon. In these 
cases, distinguishing the planet's spectrum from that of 
the star requires different techniques that don't rely on 
imaging. 

Fortunately, it has been shown recently by Charbon- 
neau et al. (2005) and Deming et al. (2005,2006) that 
the Spitzer infrared space telescope can discern changes 
in the summed light of a transiting EGP and primary 
star due to the occultation of the planet by the star dur- 
ing secondary eclipse (phase angle, a, near 0°). The 
difference in the summed light just before and during 
planetary eclipse provides a measure of the irradiated 
planet's emissions in the Spitzer IRAC bands at 3.6 /im , 
4.5 /xm, 5.8 fim, and 8.0 fim, in the MIPS band at 
24 /im , and via the Spitzer/IRS. To date, nine transit- 
ing EGPs have been discovered (Charbonneau, Brown, 
Burrows, & Laughlin 2006), four (HD209458b, TrES-1, 
HD189733b, and HD149026b) are close enough to at- 
tempt secondary eclipse measurements with adequate 
signals-to-noise, and, as of this writing, eclipses for three 
transiting EGPs have in fact been detected 2 . The corre- 
sponding planet/star flux ratios 3 at superior conjunction 

2 The planet/star flux contrast ratio for HD149026b may, how- 
ever, be too low for a successful Sptizer campaign. 

3 actually, detected electron ratios 
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are for TrES-1 0.00066±0.00013 and 0.00225±0.00036 at 
4.5 /Ltmand 8.0 ^m, respectively (Charbonneau et al. 
2005), for HD209458b 0.0026±0.00045 at 24 /zm (Deming 
et al. 2005), and for HD189733b 0.0055 ±0.00017 at -16 
/itm in the IRS peak- up band (Deming et al. 2006). Along 
with the inferences using HST/STIS of the presence of 
sodium in the atmosphere of HD209458b (Charbonneau 
et al. 2002; Fortney et al. 2003; Allard et al. 2003) 
and of photolytic atomic hydrogen in its wind (Vidal- 
Madjar et al. 2003), these data are the first direct "spec- 
tral" measurements of extrasolar planet atmospheres and 
are early harbingers of the numerous programs of EGP 
remote sensing from the ground and from space being 
planned and/or proposed. 

The secondary eclipse data for HD209458b and TrES-1 
have been subjected to preliminary theoretical analysis 
by four groups. Burrows, Hubeny, & Sudarsky (2005) 
concluded that these data are best interpreted with at- 
mospheres containing water and carbon monoxide for 
which redistribution to the night side is significant, but 
partial. They conclude that the metallicity dependence 
is very weak and predict that the flux at 3.6 /mi is higher 
than that at 4.5 /mi . They also predict a broad peak near 
10 /mi , not so obvious in the theoretical results of oth- 
ers. Seager et al. (2005) emphasize the potential effects 
of non-solar C/O ratios above 1.0, in particular the asso- 
ciated lowering of the water abundance and weakening of 
the water absorption features. They also suggest that the 
dayside reradiates most of the stellar heat absorbed and 
incorporate into their arguments the upper limit near 2.2 
/mi found for HD209458b by Richardson, Deming, & Sea- 
ger (2003). Fortney et al. (2005) have trouble fitting the 
steep spectral slope seen in TrES-1 between 4.5 /mi and 
8.0 /mi , without a significant enhancement in metallic- 
ity. With enhancements of 3 to 5, they fit the two TrES- 
1 data points to within 2-a (8 /im) and 1-cr (4.5 lira). 
Furthermore, for both TrES-1 and HD209458b they pre- 
fer uniform reradiation of the absorbed stellar light over 
the entire planetary sphere and, hence, complete heat 
redistribution. Barman et al. (2005) calculate 2D plan- 
etary atmospheres and a set of light curves for TrES-1 
and HD209458b and redistribute heat from the dayside 
with a redistribution factor / (Burrows et al. 2000; Bur- 
rows, Sudarsky, & Hubbard 2003), also used by Burrows, 
Hubeny, & Sudarsky (2005) and Fortney et al. (2005). 
They assume nightside core fluxes consistent with fixed 
values of effective temperature (T c g ) of 225 K and 500 K 
for TrES-1 and HD209458b, respectively. Barman et al. 
(2005) conclude that some redistribution must be occur- 
ring in the atmospheres of both HD209458b and TrES-1, 
but have trouble simultaneously fitting the two TrES-1 
data points. All groups find that the atmospheric tem- 
peratures are, as expected, hot and above ~1000 K, but 
the predicted planet/star contrast ratio spectra at supe- 
rior conjunction vary perceptibly from group to group. 

With the IRS peak-up measurement of HD189733b at 
~16 fira , and the secondary eclipse data anticipated in 
the near future for the remaining combinations of object 
and Spitzer band 4 , the time is ripe for a new set of the- 
oretical spectral models and predictions for HD209458b, 
TrES-1, and HD189733b at superior conjunction (a = 0), 

4 bringing the total number of data points or constraints around 
secondary eclipse to 18 (!) 



as well as for the corresponding light curves for the gen- 
eral phase angle, a. In this paper, we provide such mod- 
els and compare to the extant data to extract physical 
information about the atmospheres of these three transit- 
ing EGPs. We also make predictions for the light curves 
as a function of wavelength and the degree of redistribu- 
tion to the nightside, and explore the metallicity depen- 
dence of the secondary eclipse predictions. To calculate 
the phase light curves, we use the 2D photon transport 
code and technique described in Sudarsky et al. (2005), 
but introduce the redistribution parameter, P„, which is 
the fraction of the stellar energy intercepted by the planet 
that is redistributed to the nightside 5 . P„ = means no 
redistribution. We calculate for a given star/planet sys- 
tem and P n both the dayside and nightside atmospheric 
temperature(T)/pressure(P) profiles and the associated 
spectra, and then for a given phase angle, a, combine the 
emissions from the two hemispheres to derive the total 
planet fluxes at the Earth for 300 wavelengths logarith- 
mically spaced from 2.5 /imto 30 /imand the correspond- 
ing planet/star flux ratios. Limb darkening effects for the 
day and night sides and planetary Bond, geometric, and 
spherical albedos for the day side are automatically de- 
rived in the calculations and are not imposed artificially. 

We have opted in this paper for the P„ parametriza- 
tion, and not the / parametrization mentioned above 
and introduced by Burrows et al. (2000), because it 
is better tied to the core issue of heat redistribution, 
and because / is definitionally tied to stellar irradiation, 
which is in fact absent on the nightside. Modeling the 
heating on the nightside with a flux at the base of the 
atmosphere that accounts for the advection of heat by 
winds seemed a bit more physical than heating the night- 
side by the ensatz of external insolation. The results 
are predictions for the three transiting EGPs as a func- 
tion of wavelength or Spitzer band and six values of P n 
({0,0.1,0.2,0.3,0.4,0.5}). We thereby derive the depen- 
dence of the planet /star flux ratio spectra upon redistri- 
bution fraction and phase angle, albeit in the context of 
a simplified meteorological model. The spectral model, 
however, is state-of-the-art. In addition, we determine 
for the three close-in EGPs the approximate P„ depen- 
dence of the day and night side T/P profiles and pay spe- 
cial attention to the temporal and phase dependence of 
the flux ratios in the IRAC bands and the 24-/miMLPS 
band during a full orbit (0° < a < 180°), not just at 
superior conjunction (secondary eclipse). In this way, 
we provide a complete set of theoretical models both for 
comparison with current data and for predicting future 
measurements. Though our baseline models are for so- 
lar metallicity, we find that the metallicity dependence, 
without clouds and for solar abundance ratios, is small 
(see ©. 

2. WHY THE MID-INFRARED IS BEST AND NEW 
THEORETICAL T/P PROFILES 

It is in the near- and mid-infrared that the planet-to- 
star flux ratio is most favorable for the direct detection of 
the light of a close-in EGP. The general theory makes this 
clear (Burrows et al. 2001; Burrows et al. 2004; Sudarsky 

5 Note that with this definition, if the Bond albedo were large 
(which is the case only for cloudy models we don't discuss in this 
paper), the P n =0.5 model would result in slightly greater IR 
fluxes from the nightside than the dayside. 
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et al. 2005; Burrows 2005), but this can be most easily 
demonstrated with a graph of the contrast ratio versus 
wavelength for a sample of the closest EGPs. Figure 
Q] depicts theoretical orbital-phase averaged (Sudarsky, 
Burrows, & Hubeny 2003) flux ratios from the optical 
through 30 /im , under the assumption of complete re- 
distribution and ignoring any possible cloud effects, for 
five of the most interesting close-in EGPs: HD209458b, 
HD189733b, TrES-1, 51 Peg b, and t Boo b. The first 
three are transiting and are discussed in more detail in 
following sections. 

Figure [T] demonstrates the advantages of working in 
the mid-IR. Not only is the stellar flux lower there, but 
absorption by nascent sodium and potassium in the op- 
tical and near-IR renders measurements at shorter wave- 
lengths problematic, even for the most sensitive space 
telescopes. To date, only upper limits have been obtained 
to the planet /star contrasts in the optical and the corre- 
sponding albedos (e.g., Rowe et al. 2006). Contrasts in 
the optical of 10~ 5 to 10 -6 are expected, improving to 
only 10~ 4 to 10~ 3 in the near-IR. 

However, as Fig. [T] suggests, in the mid-IR the 
planet/star contrasts rise above 10~ 3 , making Spitzer, 
and later JWST, the preferred platforms for the spectral 
study of close-in and transiting giant planets. Note that 
there is a full order-of-magnitude range in the predicted 
average contrast ratio in the mid-IR, reflecting predomi- 
nantly the range of stellar fluxes and orbital separations. 
Fortuitously, the transiting EGPs upon which we focus 
in this paper have the most favorable ratios of the set. 

Table [T] lists the physical parameters we assume for the 
three transiting EGPs and their primaries. The most im- 
portant are the orbital radius (a), period (P), planetary 
radius (R p ), and planetary mass (M^), along with the 
stellar properties that determine the stellar luminosity 
and spectrum. In addition, Table [T] gives the surface 
gravity (g p ) employed for each EGP. The stellar irradi- 
ation spectra are taken from Kurucz (1994) and we use 
the ID and 2D atmosphere and irradiation codes outlined 
in Burrows, Sudarsky, & Hubeny (2004) and Sudarsky, 
Burrows, & Hubeny (2005). 

Figures [3j and g] depict the T/P profiles for both 
the day (solid) and night (dashed) sides of our models 
of HD209458b, TrES-1, and HD189733b for values of P„ 
from 0.0 (or 0.1) to 0.5, in steps of 0.1. P„ — implies no 
redistribution to the nightside and P„ = 0.5 assumes that 
50% of the stellar energy intercepted by the planet is ad- 
verted to the night side, where it is radiated. The thicker 
line of the two hemisphere sequences for each EGP is for 
P„ = 0.5; as Figs. [H [3J and Q] indicate, the dayside at- 
mospheres are hotter for lower P„s, while the nightside 
atmospheres are cooler for lower P„s, as would be ex- 
pected. Since for P„ = the nightside atmosphere is so 
much cooler than for the other values, P„ = nightside 
T/P profiles are not shown. 

Without clouds 6 , the albedos on the dayside are due 
solely to Rayleigh scattering. As a consequence, the 
Bond albedos are low (< 0.1; Sudarsky, Burrows, & Pinto 

6 In constructing the models present in this paper we have ne- 
glected clouds. This is sensible for TrES-1 and HD189733b, due to 
the fact that their T/P profiles do not intercept at altitude the sili- 
cate or iron condensation curves. However, for HD209458b, clouds 
in its upper atmosphere may well have some effect. We plan to 
revisit this issue in a following paper. 



2000) and P„ = 0.5 corresponds approximately to sym- 
metric losses on the day and night sides. We point out 
that the roughly isothermal region on the dayside is an 
effect of stellar irradiation, which is absent on the night- 
side. This is the region where the outward core flux and 
the inwardly penetrating stellar flux are in rough balance 
(see Hubeny, Burrows, & Sudarsky 2003 for a more de- 
tailed explanation). It is important to note that despite 
the fact that the inner boundary condition for EGP at- 
mospheres on both the day and the night sides should 
be that they have the same interior entropy and surface 
gravity, the nightside and dayside atmospheres depicted 
in Figs. [3J and H] clearly have different interior en- 
tropies. This is because the advection of heat to the 
nightside by winds and jet streams is an additional heat 
source. For the purposes of this study, we place this 
nightside source at the base of the nightside atmosphere, 
which as the Figures show is at much lower pressures 
than the convective zone on the dayside. Given that we 
have no model for the 3D general circulation, this is a 
reasonable ansatz and gives the proper emergent fluxes. 
The associated nightside T c ff is given by the formula: 

T 4 cfi = P„OT* ff (star) + T 4 cB (intrinsic) , (1) 

where f2 is the stellar flux dilution factor at the planet 
surface, T c ff (star) is the star's effective temperature, and 
cT* ff (intrinsic) is the planet's core flux, a is the Stefan- 
Boltzmann constant and T c g (intrinsic) is assumed to be 
50 K for the three EGPs. At this value of T c ff (intrinsic), 
and for reasonable values (50-75 K), its effects on the 
emergent planet fluxes and the planet/star flux ratios 
are negligible, except on the nightside for P„ = 0.0. 

As Figs. [U 02 and 2] indicate, the atmosphere 
of HD209458b is hotter than those of TrES-1 and 
HD189733b. Furthermore, the atmospheric profiles 
of TrES-1 and HD189733b are quite similar. On 
HD209458b's dayside and at a P of ~1 bar, T ranges be- 
tween 1800 K and 2100 K, while the corresponding Ts of 
TrES-1 and HD189733b range between 1500 K and 1900 
K. At pressures of ~0.1 bar, the difference in T between 
the night and day sides of HD209458b ranges from +800 
K to -200 K for P n = 0.1 and P n = 0.5, respectively, 
representing a wide, but reasonable, range (Showman k, 
Guillot 2002; Guillot & Showman 2002; Iro, Bezard, & 
Guillot 2005). The corresponding numbers for TrES-1 
and HD189733b are a bit smaller, +500 K to -200 K. 
For all three transiting EGPs, one temperature does not 
represent the atmosphere very well. In particular, the 
temperatures at the wavelength-dependent photospheres 
range by more than 500 K. Hence, one should be cautious 
when inverting a "brightness" temperature to obtain an 
"atmospheric" temperature. 

3. ORBITAL VARIATION OF PLANET/STAR FLUX 
RATIOS IN THE SPITZER BANDS 

Given these atmospheric T/P profiles, we turn next to 
the corresponding planet/star flux ratios in the Sptizer 
bands as a function of orbital phase. Figures El 
and [5] portray the light curves in the IRAC-1 (3.6 /xm), 
IRAC-2 (4.5 fim), IRAC-3 (5.8 fim), IRAC-4 (8.0 fim), 
and MIPS (24 fan) bands for HD189733b, HD209458b, 
and TrES-1. For simplicity, we have assumed orbital 
inclinations of 90° for all three EGPs and we do not 
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excise the tens of minutes during which the planet is ac- 
tually occulted and the ratio would drop to zero. To con- 
struct our models, we have used the latest EGP radii for 
HD189733b (Bakos et al. 2006) and HD209458b (Knut- 
son et al. 2006), both of which have trended down since 
discovery. For each EGP, the secondary eclipse begins 
near a ~ 7°, which due to the flatness of the curves 
near a — has the same planet/star ratios as at a = 
(superior conjunction). Figures [5] through [5] provide the 
entire orbital evolution of the contrast ratios of all three 
transiting planets, not just the values during secondary 
eclipse. The plots are versus orbit phase, or fraction of 
an orbital period, and the individual periods are given 
in Table [1] and on the figures. As expected, the high- 
est contrast ratios at secondary eclipse are for the P n 
= 0.0 models, which also manifest the largest temporal 
variation. In addition, for a given P ra , the highest values 
are for HD189733b, and in the four IRAC bands those 
for TrES-1 and HD209458b are similar. For the MIPS 
band, the TrES-1 values exceed those for HD209458b by 
—20-30%. For P n = 0.0, the variation with phase about 
the orbit mean for all five bands is about a factor of two, 
while that for P„ = 0.5 is —5-10%. There is temporal 
variation even for the P„ =0.5 models because the T/P 
profiles on the day and night sides are different, despite 
comparable total IR reradiation, and because the dayside 
optical albedos are non-zero. 

As Figs. [5] through [5] indicate, the planet/star flux ra- 
tios in the IRAC bands vary with P„ by about 30-60%, 
depending upon object and band, but only by about 20- 
30% in the 24-jumMIPS band. The MIPS band shows 
less variation with P„ at secondary eclipse because a 
change of P n from 0.0 to 0.5 represents a factor of —2 
decrease in the heating of the dayside and this translates 
into a ~ 2 1 / 4 — 1 —20% variation in the "emission" tem- 
perature. Since 24 microns is on the Rayleigh- Jeans tail, 
this is also the corresponding flux variation. 

Figure [9] shows that the predicted flux ratio in the 
MIPS band at secondary eclipse varies for HD189733b 
between —0.55% and ~0.7%, while those for TreS-1 and 
HD209458b vary from -0.4% to -0.47% and -0.3% to 
—0.4%, respectively. These large values, very different 
from the low values anticipated in the optical (Fig. [T]), 
have motivated (see also Burrows, Sudarsky, & Hubeny 
2003) the ongoing Spitzer observing campaigns. 

4. PLANET/STAR FLUX RATIO SPECTRA AS A 
FUNCTION OF PHASE AND REDISTRIBUTION 

Expanding from the narrower focus on the IRAC and 
MIPS bands in g3j Figs. EH QU and [12] portray the 
planet/star contrast spectra from 3.0 /zmto —27 /imfor 
a representative subset of P„ values and for HD209458b, 
TrES-1, and HD189733b, respectively. Shown are these 
spectra for four different phase angles (0°, 60°, 120°, and 
180°) and for the three P„s. These results are conso- 
nant with the band light curves given in Figs. O to O 
but render the full spectral dependence on P„, transit- 
ing EGP, and orbital phase. Generically, we find peaks 
near —4 /zmand —10 /imand steeper slopes from 5 /imto 
10 /zmnear superior conjunction. For P„ = 0.5, the 4- 
/impeak to 10- /im peak flux ratio increases with increas- 
ing phase angle away from a = 0, but for P n = it 
generally decreases. From 14 /imto 30 /im, the con- 
trast spectra are rather flat, implying that the 16-/imIRS 



peak-up and the 24-/imMIPS numbers should be com- 
parable. At 60° and 120°, the models with different P n s, 
while different, are most similar. Naturally, the spectra 
at 180° are entirely due to the nightside. At a given 
wavelength, the magnitudes of the phase and temporal 
variations are similar to those discussed in [J3] for those 
IRAC or MIPS bands that are closest in wavelength. Fig- 
ures [5] through [12] summarize our theoretical calculations 
and predictions. 

5. COMPARISON OF THEORY WITH DATA AT 
SECONDARY ECLIPSE AND CONCLUSIONS 

We now turn to specific comparisons between the ex- 
tant Spitzer data and our theoretical results. Figure [T3l 
portrays the planet /star flux ratios versus wavelength at 
superior conjunction for P„ = 0.5 (—complete redistribu- 
tion) and the three transiting EGPs: HD189733b (blue), 
TrES-1 (red), and HD209458b (green). Superposed as 
large squares with 1-a flux error bars are the four sec- 
ondary eclipse measurements to date (two for TrES-1, 
one for HD209458b, and one for HD189733b). Also in- 
cluded as round dots in the appropriate color are the 
band-integrated detected electron ratio predictions for 
these EGP models, with approximate band widths indi- 
cated and no error bars in the flux direction. 

A comparison between the measured points and corre- 
sponding theoretical points for this P„ = 0.5 model is en- 
couraging, particularly for the TrES-1 data at 4.5 /imand 
8.0 /im, but also for the HD189733b IRS peak-up data 
point near 16 /im. The 24-/xm point for HD209458b is 
within about 1-cr, but slightly below the theory. Since 
these P ra =0.5 models yield the lowest theoretical values 
for the contrast ratios among the set from 0.0 to 0.5, for 
the HD209458b 24- /im point this may be the best we can 
do currently. For the HD189733b point at 16 /im , the en- 
tire P„ range studied would still be consistent to within 
the 1-a range quoted (see the top left panel of Fig. [P2]) . 
with perhaps only the P„ = 0.0 model mildly discounted, 
rendering problematic for this EGP the constraint on the 
degree of heat redistribution from this one point alone. 
As Fig. [6] suggests, values of P„ of 0.0, 0.1, 0.2, and 0.3 
do not fit the TrES-1 data point at 4.5 /imto at least 2-cr, 
while values of P„ of 0.0 and 0.1 do not fit the TrES-1 
data point at 8.0 /imto the 2-cr level. The other values of 
P„ can not be excluded. Hence, we conclude that while 
some heat redistribution by winds to the nightside is def- 
initely indicated for TrES-1 and HD209458b, the degree 
of redistribution is harder to constrain, with a slight bias 
towards the larger values of P„. 

The steep slope from 4.5 /im to 8.0 /imis best explained 
by the rise to a peak near —10 /im , which in our models 
is a natural consequence of the relative strength and po- 
sitions of water bands longward of —5.5 /im(Fig. Hip . 
Note that without a large water abundance, none of 
the data nor their ratios would make collective sense. 
This was the conclusion of Burrows, Hubeny, & Su- 
darsky (2005), which we reconfirm here. Furthermore, 
the presence of the 4.67 /mi band of CO is indicated 
by the depth of the 4. 5-/im feature of TrES-1, but due 
to the fact that this IRAC-2 band measurement per- 
force sums over steeply rising fluxes in regions of the 
spectrum that bracket the 4.67 /im feature, and the fact 
that with reasonable abundances the band is saturated, 
almost nothing can be said about the CO abundance 
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(Burrows, Hubeny, & Sudarsky 2005). We predict a 
rise from IRAC-2 to IRAC-1 for all our models, indica- 
tive of the peak we generically see just shortward of 4.0 
fim . We also predict a slight peak around 10 /im , and 
a plateau from ~14 /imto 30 fim. The peak near 10 
fim might be discernible for HD189733b using the full 
capability of Spitzer/IRS. The predicted plateau seems 
suggested by the comparison between theory and data 
for the HD189733b 16-^mand HD209458b MIPS points, 
taken together, but mixing objects (as we have been 
forced to do with only four data points) is not very sat- 
isfying. 

We have calculated a P„ = 0.5 model for HD209458b 
with 10 x solar metallicity and, contrary to the conclusion 
of Fortney et al. (2005), we find that the band contrast 
ratios are within ^5% of those with solar abundances. 
This is because, without clouds, the Bond albedos are 
very low (<5%). Since changing the metallicity does not 
change the total stellar light intercepted by the planet 
for a given planet radius, the characteristic atmospheric 
temperatures are similar. What is more, we find that 
the T/P profiles are also similar, with the result that the 
fluxes and contrast ratios are little altered. We have not 
been able to trace the origin of the difference between our 
results for higher metallicities and those of Fortney et al. 
(2005). However, we interpret the very weak metallicity 
dependence of the contrast ratios at secondary eclipse for 
EGP models without clouds that we find theoretically to 
indicate that the metallicity may well be supersolar and 
large. However, by the same token, we conclude that the 
metallicity can not easily be constrained nor measured 
by secondary eclipse data alone. Cloud models, which 
we expect may be relevant for HD209458b alone among 
the three EGPs (Fortney et al. 2003; Sudarsky, Burrows, 
& Hubeny 2003), may well change this conclusion and 
variations in the C/O ratio, while we do not see any 
need at this time to invoke them to fit the four Spitzer 
data points, are still of interest (Seager et al. 2005). 

One way to significantly alter the planet /star contrast 
ratios is to introduce at altitude a strong absorber in 
the optical and near-UV, where the incident stellar flux 
can be large. In this way, the upper atmosphere is 
heated. The associated reradiated optical flux is also 
greater and the T/P profile manifests a "stratospheric" 
inversion (Hubeny, Burrows, & Sudarsky 2003). Since 
the mid-IR fluxes originate higher up in the atmosphere 
than where 7R OSSe i an d~l, the associated brightness tem- 
peratures from 4 fimto 30 /jmare also enhanced. The 
increase in the emergent fluxes in the optical and mid-IR 
leads to a corresponding suppression in the near-IR (~l-4 
fim). This potential mechanism for altering our baseline 
predictions in the Spitzer bands and for suppressing flux 
in the near-IR, particularly in the Z, J, and K bands, 
should be borne in mind. Figure [TH constructed from 
a theoretical model found in Hubeny, Burrows, & Su- 
darsky (2003), depicts an extreme version of this effect 
for P„ = models of OGLE-TR56b with (red curve) 
and without (blue curve) TiO and VO in its upper at- 
mosphere. In fact, we expect that TiO and VO are both 
flushed out of the upper atmosphere by the coldtrap ef- 
fect, but suppressing this effect allows us to make the 
general point. Note that the bumps near 10 /imand 4 
fimseen. in our fiducial models (Figs. [TOl - [T2l and Fig. 
[T3|) can be altered, shifted (4 fim), or muted (10 fim) by 



this upper-atmosphere absorption effect, so if we fail to 
see these features as predicted interesting stratospheric 
or upper-atmosphere chemistry might be implied. Con- 
versely, the presence of these bumps will put useful limits 
on such upper-atmosphere absorbers. 

As mentioned in SJU we employ the latest measure- 
ments of the transit radii of each of the three EGPs in de- 
termining the planet/star contrast levels. However, these 
radii, being transit radii that probe along the chord of 
the planet in the optical (Burrows, Sudarsky, & Hubbard 
2003), are not strictly the appropriate radii to use in de- 
termining the planet/star contrast ratios. They are close, 
but the total stellar energy intercepted by the planet de- 
pends upon wavelength and is different in the near-IR 
water bands, where the transit radius should be slightly 
larger than in the optical (for HD209458b, by ~2-3%; 
Fortney et al. 2003). Furthermore, planetary emission is 
from a radius that is not corrected for by the "transit ra- 
dius effect." For HD20948b, this radius difference can be 
8-10% (Burrows, Sudarsky, & Hubbard 2003), while for 
TrES-1 and HD189733b, due to the lower atmospheric 
temperatures and higher gravities the effect is smaller 
(<5%). The upshot is that the predicted flux ratios for 
HD209458b could be smaller by as much as ~15%, bring- 
ing the 24-/imMIPS point better in line with our predic- 
tion, but introducing further ambiguities into predictions 
at all wavelengths until the radius issue is resolved. No 
one doing theoretical secondary eclipse calculations has 
yet corrected for these subtle radius effects. Moreover, 
even the measured transit radii retain a residual ambi- 
guity due to the systematic uncertainty in the stellar ra- 
dius, which could easily be 5-10%. Knutson et al. (2006) 
argue that measurements of the star, constraints of stel- 
lar evolution theory (Cody & Sasselov 2002), and the 
detailed fits to the HST/STIS transit measurements to- 
gether yield a transit radius for HD209458b with an error 
of only ~2%. Perhaps, but our same concerns apply to 
TrES-1 and HD189733b. In sum, ambiguities in the ap- 
propriate radii to employ in comparing secondary eclipse 
data with theory remain and slightly compromise their 
interpretation. 

However, that the data and theory we have devel- 
oped here correspond as well as they do is gratifying. 
In fact, the theory does a good job fitting the four 
secondary eclipse data points (Fig. [T2]), whatever the 
ambiguities. In addition, there is evidence for redis- 
tribtuion to the nightside, particularly for TrES-1 and 
HD209458b, (though its specific magnitude remains to 
be determined), and the presence of H2O and of CO is 
strongly indicated. Moreover, we find that the metallic- 
ity dependence of cloud-free models is quite mild, but 
that ambiguities in the radii remain to slightly compro- 
mise the interpretation of the data. Due to the greater 
sensitivity in the IRAC bands to variations in P„ (Figs. 
[5]-|8]), such data have greater potential to determine the 
degree(s) of redistribution. Data off secondary eclipse at 
other phase angles, and particularly in IRAC-4 (Figs. ITOl 
- f]~2]) - would further constrain the models, but the most 
propitious phase angles in this regard are larger than 90° 
and, hence, will prove very difficult to measure. How- 
ever, JWST, with its two to three orders-of-magnitude 
greater sensitivity in the mid-IR, will be able to mea- 
sure a large fraction of the planetary light curves. In the 
shorter term, data, or even upper limits, at the 14 other 
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anticipated Spitzer band points are anxiously awaited. 

We thank Christopher Sharp, Bill Hubbard, Jaymie 
Matthews, Jonathan Fortney, Mike Cushing, and Drew 
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available in electronic form from the first author upon 
request. A web-based calculator for determining the op- 
tical and near-IR light curves of EGPs in wide orbits as 
a function of orbital phase, wavelength, semi-major axis, 

orbital inclination, and eccentri city is now av ailable at 

|http:/ /zenith.as.arizona.edu/~burrows/phase/lightcur ve.php[ 
A similar calculator for close-in EGPs is under develop- 
ment. 
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TABLE 1 

Reference Data for Transiting EGPs 



EGP Star M* (M Q ) R, (R ) a (AU) P (days) M p (Mj ) R p (Rj) Log 10 g p (cgs) [Fe/H] 

HD209458b GOV 1.05 1.1 0.045 3.524 0.69 1.32 2.99 0.01 

HD189733b K2V 0.82 0.75 0.031 2.22 1.15 1.154 3.25 -0.03 

TrES-1 KOV 0.87 0.86 0.0393 3.03 0.76 1.08 3.21 0.0 



References. - Rj (Jupiter's r adius) — 7.149 x 10 4 km; see J. Schneider's E xtrasolar Planet Encyclopaedia at 

http:/ /www. obspm.fr/cncycl/cncycl. html and the Carnegie/California compilation at http://exoplanets.org 
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Fig. 1. — Logarithm of the planet/star flux density contrast ratio versus wavelength from 0.5 /imto 30 /an, for five close-in EGPs: 
HD189733b, HD209458b, TrES-1, 51 Peg b, and r Boo b. These ratios have been averaged over a full orbit, assume complete redistribution 
of heat using the old prescription (/ = 0.25; Burrows, Sudarsky, & Hubbard 2003), and are from cloud-free models. Note the variation of 
more than three orders-of-magnitude from the optical to the mid-IR, where, all else being equal, the ratios are most favorable and Spitzer 
is currently the platform of choice. See text for a discussion. 



9 




1000 2000 3000 



T [K] 



Fig. 2. — Theoretical profiles of the logarithm base ten of the pressure (P, in bars) versus temperature (T, in Kelvin) for both the daysidc 
(solid) and the nightsidc (dashed) atmospheric models of HD209458b, for P n s from 0.0 to 0.5 in steps of 0.1. The P„ = model for the 
nightsidc is omitted. The thicker lines are for the P n = 0.5 models. See text for details and discussion. 



10 




1000 2000 3000 



T [K] 



Fig. 3.— Same as Fig. [2] but for TrES-1. 
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Fig. 4.— Same as Figs. |2]and[3] but for HD189733b. 
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Fig. 5. — Light curves versus orbital phase (or orbit fraction) of the planet/star flux ratios in the IRAC-1 band centered near 3.6 /im, 
for the transiting EGPs HD189733b (blue), HD209458b (green), and TrES-1 (red) and for P n s from 0.1 to 0.5, in steps of 0.1. For a given 
object or color, the higher curves are for the lower values of P n and the dependence upon P n is monotonic. Indicated by vertical dashed 
lines are the positions of the a = 0° (superior conjunction, full phase) and a = 180° phases. Given in the legend are the orbital periods of 
the three EGPs. See text for a more detailed discussion. 
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Fig. 6. — Same as Fig. \5\ but for the IRAC-2 band centered near 4.5 fim . 
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Fig. 7. — Same as Fig. \5\ but for the IRAC-3 band centered near 5.8 fim. 
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Same as Fig. \5\ but for the IRAC-4 band centered near 8.0 (im 
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Fig. 9. — Same as Fig. [5] but for the MIPS band centered near 24 /im. 
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Fig. 10. — The HD209458 planet/star flux ratio spectra for three representative P n s (0.1, 0.3, 0.5) and four phases (a = 
0°, 60°, 120°, andl80°). The wavelength range is 3 /imto 27 /an. See text for a short discussion of the various features and system- 
atics and Figs. [5] through [9] for related theoretical results focussed on the IRAC and MIPS bands. 
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Fig. 11 — 



Same as Fig. 1101 but for TrES-1 and its primary. 
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Fig. 12 — 



Same as Fig. \W\ but for the HD189733 EGP/star system. 
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Fig. 13. — Theoretical planet/star flux ratio spectra during secondary eclipse for wavelengths from 3 fim to 27 fim , for P„ = 0.5 (^complete 
heat redistribution) and for the three EGPs HD189733b (blue), TrES-1 (red), and HD209458b (green). Values for 300 wavelength points 
spaced logarithmically are included in each curve. Superposed as large squares with 1-a error bars and in the appropriate color are the 
current measurements at secondary eclipse (Deming et al. 2005, 2006; Charbonneau et al. 2005). The smaller filled circles are the 
predictions for the four IRAC bands, the MIPS band at 24-/im , and the IRS peak-up band near 16 (im . To derive these numbers we 
have performed a band-average of the ratio of the detected electrons. One should compare the small dots with the large squares of the 
same color and in the same wavelength band to draw conclusions about the fidelity of the models and the character of the associated EGP 
atmospheres. See the text in j5]for a discussion of some of the inferences drawn from the current dataset. 
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Fig. 14. — Two examples of theoretical flux spectra (F„, in erg cm" 2 s" 1 Hz -1 ) from the surface of the close-in EGP OGLE-TR56b 
from the optical to 30 /im, with (red line) and without (blue line) TiO and VO in its upper atmosphere. For both models, P n has been set 
equal to 0.0. The differences due to the possible effects of a strong optical absorber placed in the upper atmosphere of a close-in, hot EGP 
in the optical, near-IR, and mid-IR are clearly seen. In the toy model with "TiO/VO" (red line), a stratospheric inversion is produced, the 
mid-IR and optical are enhanced, and the near-IR is suppressed. Model numbers were taken from Hubcny, Burrows, & Sudarsky (2003). 
See text for a discussion. 



